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ABSTRACT 

It is shown that the results of maximum likelihood treatment of Monte 
Carlo simulation with constant production rate of 7*6 SNU and 1*8 SNU 
are consistent with the constant production rate when the tests of hypo- 
theses(e*g* t-te st ,3C -te at, Wilcoxon-Mann- Whitney test, run test etc*)has 
been applied to the two groups of data formed from sunspot minimu m range 
and sunspot maximum range whereas the real data pulsates with the solar 
activity cycle*It is shown that SN flux-change is in opposite phase to 
the solar activity cycle and lags behind the latter ty about one year* A 
correlation between SN flux and the cosmic rays has also been suggested* 

Reychaudhuri(1,2,3,4*5,6)and Gavrin et al( 7) showed that the solar neu- 
trino( SN^flux data is varying with the solar activity cycle • Ray chaudhur i 
(3,4, 5, 6) showed that the fluxes will be higher in the sunspot minimum 
range than in the sihksppttmaximum range* This is supported ty the data 
presented by Rowley et al(8)at the solar neutrino and neutrino astronomy 
conference held at Lead, South Dakota* Apart from the above variation Ray- 
chaudhuri (5,6) showed that there is also a quasibienniel variation from 
1970-1975 and 1979- 1984* Ibis type of variation is already observed ty 
Ctoevyshev(9)and fllisetti( 10) in the cosmic ray intensity (proton flare)* 

In this paper an estimate of the moving average(5 successive run numbej) 
data is considered to find the statistical significant variation of SN 
flux with the solar activity cycle of 11 years by t-te at pc*— test, wucoxon 
•* Mann- Whitney test and run test* To search variation in a convincing man- 
ner the above tests of hypotheses have been applied to the data that had 
been generated by Monte Carlo simulation and background parameters are 
typical of those in the actual experiment* 

SOLAR ACTIVITY CYCLE AKD(a ^MOVING AVERAGE DATA. ( b )M0NTB CARLO SIMULATED 
DATA WITH CONSTACT PRODIKTION RATE( ifa ^ Sm? AND(iiT^T' SNlT “ 

Here we apply the same proceedure as in Raychaudhuri( 5, 6 )• 
l)Student t-teat;-a(i) Moving average data 

We form the two groups,the first group can prises sunspot minimum range 
(about 4*6 years)and the second group i*e*,the rest c emprises sunspot 
maximum range* Let us now cdlect the data from the run number 36 to 58 
from August 1974 to February 1979(about 4*6 years)for the first group* If 
we now set up the null hypothesis H : m, « m a against the alternatives m,> 
m-t *The statistics t i s gi ven by (he re the difference of variance is not 
significant)*. , _ pr- Jn/Wx - 

r - >lv >/ N(+w | 4 \p^s*+7i7s?- 

where * N, -fN^-J^In the above case N, =23,N i =27, < > ) =48 which gives t=4*54, 
Thus we can conclude that 99*9# of the data pulsates with the solar acti 
vity cycle* The standard error of the difference of means is 4*54€* H ence 
the data are with the asswption that the two means are 

eoual* 

(ii) Yearly average of the moving average dataj-Here one group comprises 
the yearly average data from 1975 to 1978 and the other grbup comprises 
. rest of the moving average data* Vb get for 9 d*f* t =2* 76 ‘The sta nda rd 
error of the difference of the means is 3* 1C which suggests that the 
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,Pig*l3?Ar production rates 
{(observed in the Chlorine SH 
^ experiment 1970-1982 and the 
Ijnoving average data (dashed 
fjcurve). 

I'Fig* 2 and 3 ‘Results of maxi- 
|mun likelihood treatment of 
Monte Carlo simulations with 

5 constant production rate of 
7*6 SHU and 1*8 SHU respect- 

— — 

means of the two samples are 

'not equal* 

b)Monte Carlo simulated data 
:-i)7*6 SHU andii)1*8 SHIM* 
form the group* as in moving 
(average data* In the case of 
7*6 SHU the difference of 
variance is significant* we 
have to use the formula for 
fro the difference of two 
means* x— y~ 

4 " 

glror d*f* 47,te1*19*Tne stand- 
ard error of the difference 
ijof two means is 1*1 6 «B»nce 
the two means are equal i*e» , 
the [steady flux is accepta- 
ble* In the case of 1*8 SHU 
(the difference of variance 
is not significant* Hence t= 
1*84,d*f* 45* The standard er- 
ror is found to be 1*8 • 

Hence the data are consist- 
ent with the assumption that 
? the steady flux is acceptab- 
“(jle^Hext we take the yearly 
average of the Monte Carlo 
simulated data from fig 2 & 

3 and form the two groups «s 
in t-test* Hare in both cases 
difference of Variance is 
not significant* In case 7*6 
SHn,t«0*43 d*f* 9* In case 

, fi1*8 SNU,t*1*01 d*f* 9*HBre 

the standard error is found to be (i)0*53 6 anU(ii) 1 € xuspectively* Hen- 
ce the data are consistent with the assumption that the steady fltifcSJjjg 
acceptable* . J 

II) , Xv -testj-(i)W0 form the group as in the t-test* R»r moving average 
data £T % =0*02 and for d*f*49,£*- =73*19* r flms the hypothesis of stdfcifrtf 
flux is not acceptable* Per the case of 7*6 SHB, (TV)* 291, X v »47* 84 d*f*47* 
again for 1*8 SHU, & =0*10, 5^46 *83 d*f* 45* Hence the difference of two 
samples are not significant for the case of 7*6 SNU and 1*8 SHU* 
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(ii)if we take the yearly average of the moving average data and Kent* 
Carlo simulated data* We fora the groups as in t-test «We get <7=*01, 5t x « 
I8*,55,d*f* 10 for moving averse data* In this case steaejy SN flux is not 
acceptabOLe*Again (i) (5' Z at*09> =9*39, d*f* 9 for 7*6 SNU and (ii)<J" , ®*C5, 

^ l s*1Q*/$,d*f* 9 for 1*6 SNIMfence the steady SN flux is acceptable ia 
Monte Carlo simulated data* 

( III )WULcoxo»-Manx>- Whitney test: -Here we take the yearly average data of 
moving average data and Monte Carlo simulated data 7 *6 SNU and 1*8 SNU 
and form the two groups as in t-test and ^-test* If we set up the m£Ll; hy- 
pothesis H 0 jthe distribution of x and y are identical*If a two sided test 
at 5J6 is desired for m=4,n»7 we use 30.425 as the accepted region*For 
this test oL =p Ho (U 6 3 ©r U >25)=*042*. 

s)Moving average data:- we write the two sets of data as y: # 3552,*4740, 

• 3546, *4007, *3089 , • 2469, • 4707 and. xj* 5305,* 4437, *7482,* 5421* Putting these 
numbers in order and keep track of which numbers came from x and which 
from y:yyyyyacyypaa*There are .2 inversions of the y 1 s and the rank of yfs 
is 1+2+3+4+5t %8 =30 =2+24£ Since the number of inversion of the y»s is 
not within 3 and 2$ ,we reject* the null hypothesis,!*®* the constant SN 
flux is not acceptable* 

b)(i)7*6 SNU :-y: 1*408,0*993, 1*430, 1*286, 1 *277,1 *928 and x: 1600,1*396, 
1*320, 1*098* Putting these numbers in order and keep track of which num- 
bers came from x and which came from y: yxyyxxyyxyy • There are 16 inversi- 
ons of the y<s and the rank of y's is 1+3+4+7>*0+1OVl1 «44» 

(ii)1*S SNU:-y : 1* 545, *600 1 , • 1 65 , *805,* 537, *307, • 378 and xs* 3755,* 4900, x* 
•j 518 ,*^581* There are 8 inversion of the y*s and the rank of y's is 36* 

Since the number of inversion of the y is within 3 and 25,we accept the 
null hypothesis i*e*, steady SN flux is acceptable in both the Monte Carlo 


cases* 

(IV) Run testj-In the case of run up and down the total number of runs is 
approximately normally distributed under the assumption of randcmness 
with mean and variance is given l y E(r) =l/3(2n-l) and V(r)«y(90)(6n-29)* 
For moving average data :— n»50,2fe)f33,V(r)*!8*57*In this case there are 
21 runs either positive or negative* If we take the critical region to be 
=o. 20 vouLd yield a K of 33-0*64 v/ 8* 57 =30* I 543*Thu3 we cannot accept 
the randomness since 21 <’30* 543* Thus we can suggest that the Obnetaat SN 
flux is not acoeptahle*For case 7*6 SNU,n»49,E(r)-32.33 and V(r)=8*389*a 7 i 
this case there are 31 runs either positive or negative* For the cage 1*8 
SNU,n=47,E(r)dO*31 and V(r)=8*03*In this caee there are 29 runs is either 
positive or negative* If ew take the critical regio n to be r<£K,a choice. 
of d =0*20 would yield a K of (1)32.33-0*84^8*389 =29*9|(ii)31-0*84V 803 
=28 *,62* Thus we accept the randomness in both the case since (i)3l > 29*9 
and (ii)29)> 28*62 *So we can conclude that the Monte Carlo simulated data 
da not pul sate with the solar activity cycle* 

Sun-root cvele and SN flux:- We have taken the yearly average smooth suns- 
Jot frem solar geophysical data,yearly average of the moving average data 
frem fig*|1*and yearly aveaage of the cosmic ray data I wxth E =0* 1-5*8 
Gev(ll)*W© have displayed the data in Table I* ' . 

The correlation coefficient between sunspot cycle and 9 Ar production 
data ar© r(S,Q M )« -0* u 6l,when^T®-1,r(S,Q rt )=-0^9 

0*35 whenAT®-2 where AT * T s - To (in years )*ttus the 3? Ar production 
rate anticorrelates with the sunspot cycle* The correlation coefficient 
between the galactic cosmic rcy intensity and 3f Ar production data r(I 0 , 
q m )=0*58 when4T=0,r(I 0 ,Q M )=0*31 when/1 T=-1,r(Ip,Q N )*0*39 when/ 
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Table 1* Yearly averages of smooth sunspot number(S), Ar production data 
( Q m ),yearly average of cosmic ray data I 0 with E=0* 1-5*8 Gev* 

1971 1972 1973 1974 1975 1976 1997 1978 1979 1980 1981 

sSyt'^bfr 66-6 68-9 38-5 34-4 15-5 12-6 27-5 92-7 155-3 155 140 

oftX ‘2ST .06 -44 -35 -40 - 53 -.44 - 75 - 54 ,31 -25 -47 

Ar production 
data( ) 

FiUottc'XaJc 8,01 2,66 2-80 2-7.1 >11 >20 3-29 2-79 1-89 1-24 1-13 
ray intensity 


4T=1» Where d Tb Tr 0 -.Ta rt (in years)* Thus there is an important correlation 
between the SN flux in Davis* data and galactic cosmic rays* This shows 
that the galactic cosmic rays and SN flux undergo the essential changes 
in the solar acivity cycle from 1970-1982* 

Erom the SN flux data we see that SN flux is higher during the ascend* 
ing phase( about 3*3 years )and descending phase( about 2-3 years )of the 
solar cycle avoiding the sunspot maximun time* We call the time 2-3 years 
before the sunspot maxxnnm as the second m-i nimup time* Thus SN flux is 1 
higher in both second minimum time and second maximum time and SN flux is 
lower in both first minimumand first sunspot maximum stage* In addition to 
two prominent maximum in the year 1957 and i960 they(9,10)herye found thi- 
rd maximum of smaller amplitude in the cosmiv ray intensity out veru clear 
in 1963 *Erom fig*.1 we see that run nunber 37 is connected with third max- 
imum* This may have some connection with the quasibienniel variation of SN 
flux found from 1970-1975 and 1979-1984* the moving average data in 
fig* 1 it appears that SN flux varies with a period of about 2*4 to 2*5 
years from 19 X- 1975* 

In conclusion it may be remarked that the present SN data exhibits a 
high level of statistically significant variation of SN flux with the 
solar activity cycle* The above results strongly suggests that the solar 
activity cycle ii due to the pulsating character of the nuclear energy 
generation inside the core of the sin* We suggest that the chlorine SN 
experiment should be continued to get the data like the sunspot data etc* 
and this data could be of great importance to our understanding of the 
sun and cosmic rays* 

Acknowledgement:- The author would like to thank Professor R* Davis,Jr* for 
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